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A STUDY OF THE EFFECT OF SELECTED MATERIAL PROPERTIES 
ON THE ABLATION PERFORMANCE OF ARTIFICIAL GRAPHITE 

By Howard G. Maahs 
Langley Research Center 

SUMMARY 

A total of 18 material properties have been measured on 45 different, commercially 
available, artificial graphites. A correlation based on 10 of these graphites was devel- 
oped, permitting the difficult -to -measure property, thermal conductivity, to be predicted 
in terms of the easier -to -measure property, electrical conductivity. Other correlations 
between pairs of the 18 graphite material properties have been developed where possible, 
reducing the original 18 material properties to 10 mutually independent properties. Abla- 
tion performance (as determined by total length change during test, recession rate, sur- 
face temperature, surface texture, and degree of gouging) was measured on the same 
45 graphites in a nominal Mach 2 airstream having a nominal enthalpy of 2.17 MJ/kg and 
a stagnation pressure of 5.6 atm. Nominal heating rates were 513 W/cm2. 

Multiple regression equations were formulated relating ablation performance to the 
reduced set of 10 independent material properties, thereby identifying those material 
properties having the most pronounced effect on ablation performance in the present test 
environment. These properties are maximum grain size, density, ash content, thermal 
conductivity, and mean pore radius. For optimization of ablation performance, grain 
size should be small, ash content low, density and thermal conductivity high, and mean 
pore radius large. Maximum grain size should be kept below 0.22 mm to insure a smooth 
graphite surface during ablation. 

The method of stepwise multiple regression employed in this study for determining 
which properties are important in affecting ablation performance has been shown to be 
both powerful and efficient. Its soundness is demonstrated by the fact that it properly 
predicts ablation performance to be related to certain material properties known from 
theoretical or experimental considerations to affect ablation performance. In order to 
determine other material properties which might also have an important effect on abla- 
tion performance, a more extensive list of properties than the 18 considered here could 
be employed. 

It is of some additional interest that, in the course of the present investigation, 
Graph-i-tite "G" graphite was found to be more erosion resistant than the usual standard 



of comparison, ATJ graphite. Of the two forms of Graph-i-tite "G," the large -grained 
form eroded with a rough surface while the fine-grained form eroded with a smooth sur- 
face. Accordingly, the fine-grained form of Graph-i-tite "G" appears attractive for 
future study. 


INTRODUCTION 

Carbon is a unique material existing in a wide variety of forms: natural graphite, 
carbon black, industrial artificial graphite, and more recently, the pyrolytic carbons and 
graphites and the glass -like carbons. Although all except natural graphite are artificial 
in the sense that they are manufactured, the term "artificial graphite" is commonly 
reserved for the conventional industrial graphite fabricated from carbonaceous filler and 
binder materials and subsequently graphitized. 

Artificial graphite has been of particular interest for some time as an ablation 
material because of its many desirable structural and thermal properties at high tempera- 
tures, as well as its machinability and ease of fabrication into relatively large pieces. 
However, a detailed understanding of this material and its performance in various appli- 
cations is severely complicated because it is fabricated in many forms by a wide variety 
of proprietary formulations and techniques. Furthermore, each manufacturer usually 
produces a broad line of artificial graphites, each intended for different applications. 
Because of these complications, there exists to date no simple, meaningful way to select, 
from the hundreds of artificial graphites manufactured, the ones best suited for a partic- 
ular use, such as for ablation materials. Also, because of the proprietary nature of these 
artificial graphites, there is even no way by which any two graphites can be determined 
to be fully equivalent, whether produced by the same manufacturer or not. 

In an attempt to resolve this problem, it was reasoned that if the primary funda- 
mental physical, chemical, and structural properties, as well as ablation performance, 
were determined for a large number of different graphites then possibly some dependence 
of ablation performance on the fundamental properties could be found. This information 
would aid in selecting particular graphites for consideration in ablation applications, but 
would be of greatest value in establishing the prime collection of properties desirable to 
fabricate into a graphite in order to optimize its ablation performance. 

Some of this basic property information is already available. Chemical impurity 
data are available in reference 1 on 40 selected grades of commercially available artifi- 
cial graphite, with the impurities grouped by reported effect on oxidation rate, that is, 
whether catalytically accelerating or inhibiting. (Ref. 1 indicates that the combined effect 
of the natural impurities in graphite should be accelerating, and this is confirmed experi- 
mentally in ref. 2 for four graphites differing greatly in impurity content, grain size, and 
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density.) Crystallographic data on crystallite sizes and lattice dimensions are presented 
in reference 3 for the same 40 artificial graphites studied in reference 1. 

The purpose of the present investigation is to measure for a wide variety of arti- 
ficial graphites those additional material properties which are suspected of having poten- 
tially large effects on ablation performance, to determine the ablation performance of 
these same graphites, and to relate ablation performance to the measured properties. 

A second purpose is to develop correlations, where possible, between pairs of the graphite 
material properties themselves. 

In this report are presented measurements of physical, chemical, and structural 
properties of 45 commercially available artificial graphites. These properties include 
density, maximum grain size, ash content, thermal conductivity, interlayer spacing of the 
filler carbon, and measurements of internal volume, surface area, and pore radii. Cor- 
relations between pairs of these properties have been made where appropriate. Also 
presented for these same graphites are the following ablation -performance data: total 
length change, stagnation -point recession rate, steady-state stagnation -point surface 
temperature, surface texture, and degree of gouging. These ablation data are analyzed 
in light of the material -property data, with multiple regression equations being formu- 
lated which relate ablation performance to them. 

SYMBOLS 


Ap pore area, m2/g 

A ash content, ppm 

a lattice dimension in the a direction, A 

ai coefficient 

bj coefficient 

c constant 

dg interlayer spacing of the filler carbon, A 

Gjn maximum grain size, mm 

i index 
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thermal conductivity, W/cm-K 


pressure, N/m2 
pore radius, /xm 
most probable pore radius, pm 

average absolute percent residual (or deviation) about regression 

pore radius in the micropore region, A 

average pore radius, A 

mean pore radius, A 

most probable micropore radius, A 

surface area, m2/g 

linear recession rate at the stagnation point, cm/sec 
stagnation -point surface temperature at steady state, K 
total pore volume, cm^/g 

cumulative volume in pores up to pore radius r, cm^/g 
pore volume, cm3/g 
i*^ material property 

a general ablation performance parameter 
degree of gouging 

total specimen length change resulting from a 30 -second exposure to the 
test environment, cm 



e 

®air 

%g 

^air 

Pp 

Prp 


total void porosity, percent 

contact angle between mercury and graphite, degrees 

air -open void porosity, percent 

mercury -open void porosity, percent 

air displacement density, g/cm^ 

bulk density, g/cm^ 

theoretical density, g/cm^ 

surface tension, N/m 

electrical conductivity, (f2-cm)"^ 

surface texture 


DESCRIPTION OF MATERIALS 

The particular artificial graphites considered in the present study were intention- 
ally selected to include a wide range of material properties. These graphites are listed 
in table I. All are commercially available. Several were obtained from the manufac- 
turers in a specially purified form (presumably purified by high-temperature halogena- 
tion) in addition to the conventionally supplied form and are noted in table I by the letters 
"GP" appended to the standard grade designation. With the exception of graphites AHDG, 
AGSX, AGOT, and the Graph-i-tite grades, which were fabricated by extrusion, all of the 
graphites were fabricated by molding. 

Specimens for the mate rial -property and ablation -performance measurements were 
drawn at random from as -supplied graphite billets obtained directly from the manufac- 
turers. The required specimens for each particular graphite were drawn from the same 
production batch, and, in all but a few cases, were drawn from a single billet. Prior to 
cutting the specimens from a billet, the billet was trimmed to eliminate the possibility of 
nonrepresentative skin effects confusing the subsequent material -property or ablation- 
performance measurements. Replicate measurements of the material properties on dif- 
ferent specimens of the same graphite yielded no unreasonable variations in measured 
values. 
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MATERIAL -PROPERTY MEASUREMENTS 


Bulk Density 

Bulk density pp was determined from the mass of a graphite specimen divided by 
its volume. Specimen mass was measured with an analytical balance, and specimen vol- 
ume was calculated from the dimensions of an accurately machined specimen. Specimen 
volume was on the order of 3 cm^. Measurements were made on three specimens of 
each graphite and the results averaged. 

Air Displacement Density 

Air displacement density p„ 2 „ was determined as the mass of a specimen divided 
by the volume of the solid in the specimen displacing air at 2 atm. Mass was measured 
by an analytical balance, and the solid volume was measured with a gas pycnometer 
(Beckman Instruments, Inc., Model 930). Three specimens of each graphite were mea- 
sured and the results averaged. 


Theoretical Density 

The theoretical density prj, is defined as the true density of the graphite crystal, 
and therefore depends on the lattice dimensions a and dg. Nightingale (ref. 4) gives 
an expression for the theoretical density p^ = 91.93/^a^dc. The lattice dimension a 
was shown in reference 3 to be a constant at 2.4614 A. Substitution of this value reduces 
the expression for theoretical density to p,p = 7.587/dc. The interlayer spacing dg 
was obtained from reference 3. 


Total Pore Volume 

The total pore volume Vp is a specific pore volume defined as the volume of 
mercury which can be intruded into the pores of the graphite specimen at a pressure 
of 103.4 MN/m2 (15 000 psia) divided by the mass of the specimen. All mercury intru- 
sion measurements were made on a mercury porosimeter and were performed by the 
American Instrument Co., Inc., on samples furnished by the National Aeronautics and 
Space Administration. Total pore volume was measured on two specimens of each 
graphite, and the results were averaged. 

Most Probable Pore Radius 

The most probable pore radius Rp is defined as that pore radius at which occurs 
the greatest volume increment of mercury intruded into the graphite per increment in 
pore radius. The experimental data are obtained with a mercury porosimeter, with the 
pore radius being related to the applied intrusion pressure P by the equation 
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R = -2a cos 0/P where a is the surface tension of mercury and 0 is the contact angle 
between mercury and graphite. Taking a= 480 x 10"^ N/m and 0 = 140° (ref. 5) 
there results R = 735.4 x lO^/P where P is in N/m^ and R is in /im. Plots of 
mercury intrusion volume as a function of intrusion pressure were prepared by the 
American Instrument Co., Inc., on samples furnished by NASA. Most probable pore radii 
Rp were obtained from these plots by the author. Since intrusion pressure ranged from 
about 13.8 kN/m2 (2 psia) to 103.4 MN/m2 (15 000 psia), the distribution of radii observ- 
able lies between approximately 0.007 jum and 53 /im. Reported results are the average 
of measurements on two samples. 


Maximum Grain Size 

The maximum grain size Gm is the maximum particle size of the filler carbon 
used in the manufacture of the artificial graphite, as conventionally determined by screen- 
ing. The maximum grain size was obtained from the manufacturers. 

Pore Volume 

Pore volume Vp is a specific volume defined as the void volume in pores having 
radii less than or equal to 300 A per unit mass of the specimen. It was determined from 
the nitrogen adsorption isotherms of the graphites by the method of Barrett, Joyner, and 
Halenda (ref. 6). Although nitrogen has a molecular diameter on the order of 3 A (ref. 7), 
as a practical matter, the lower limit of the pore radii included in the volume measure- 
ment by this method is about 7 A. The generation of the adsorption isotherms and the 
calculation of pore volume were performed by the American Instrument Co., Inc., on 
specimens supplied by NASA. Measurements for each graphite were made on two speci- 
mens and the results averaged. 


Pore Area 

The pore area Ap is a specific surface area defined as the surface area in pores 
having radii less than or equal to 300 A per unit mass of the specimen. Pore area was 
calculated from the same nitrogen adsorption isotherms used to determine pore volume 
Vp by the American Instrument Co., Inc., according to the method of Barrett, Joyner, 
and Halenda (ref. 6). Measurements were made on two specimens and the results 
averaged. 


Surface Area 

Surface area S is the total solid surface area per unit mass as calculated from 
the first part of the nitrogen adsorption isotherm by the classical BET method (ref. 4) . 
The same nitrogen adsorption isotherms used to determine vp and Ap were used and 
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the American Instrument Co., Inc., performed the calculations. Results are the average 
of individual determinations on two separate specimens. 

Average Pore Radius 

The average pore radius ra is defined as that pore radius which divides the pore 
volume distribution into two equal parts (ref. 8). It was calculated from the nitrogen 
adsorption isotherms of the specimens by the American Instrument Co., Inc. Measure- 
ments were made on two separate specimens and the results averaged. 

Most Probable Micropore Radius 

The most probable micropore radius rp is defined as that pore radius at which 
occurs the greatest change of cumulative pore volume v per change in pore radius r. 
It is analogous to the most probable pore radius Rp with the difference being that the 
most probable micropore radius is determined from the nitrogen adsorption isotherms, 
whereas the most probable pore radius is determined from mercury porosimetry. 

(Recall that mercury porosimetry yields a pore radius distribution between 0.007 jim and 
53 lira, and that nitrogen adsorption yields a radius distribution between 7 A (0.0007 /im) 
and 300 A (0.030 lira).) The most probable micropore radius was determined by the 
author from plots of dv/dr as a function of r prepared from the nitrogen adsorption 
isotherms of the specimens generated by the American Instrument Co., Inc. Results are 
the average of individual determinations on two separate specimens. 

Mean Pore Radius 

The mean pore radius rj^, an alternate definition of a characteristic pore dimen- 
sion, is defined as 2vpy/S, where Vp is the pore volume, and S is the surface area 
(ref. 8). It was computed from Vp and S as determined from the nitrogen adsorption 
isotherms. The results represent an average of two measurements, one on each of two 
samples. 


Total Void Porosity 

The total void porosity e is defined as the ratio of total void volume — whether 
open (i.e., accessible to air) or closed (i.e., inaccessible to air) - to bulk volume. It is 
expressible in terms of the densities already discussed and was calculated from the for- 
mula e = 1 - Ppy/Prp. 


Air-Open Void Porosity 


Air -open void porosity 
bulk volume. It was calculated from 6. 


^air defined as the ratio of void volume open to air to 


air 


= 1 - 


Pp/Pai 


air 
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Mercury -Open Void Porosity 


Mercury -open void porosity alternate definition of open void porosity, is 

defined as the ratio of the void volume as measured by mercury intrusion to bulk volume. 
It differs from air -open void porosity in that more pore volume should be open to air than 
to mercury (in spite of the high mercury intrusion pressures). It was calculated from 
the mercury porosimetry data by the formula 0jjg = PpVp. 

Ash Content 

Ash content A is a relative measure of the concentration of chemical impurities 
present. Experimentally, it is the mass of residue remaining after complete combustion 
of the graphite, expressed in mass parts per million of the original graphite. Ash con- 
tent is a convenient measure of the relative concentration of total elemental impurities 
in an artificial graphite (ref. 1). Except where noted, ash content was obtained from 
reference 1 where a breakdown of the total chemical impurities in the graphites into 
separate chemical elements is also given. 

Interlayer Spacing 

Interlayer spacing of the filler carbon dc is a measure of the spacing between 
layer planes in the filler -carbon crystal lattice, and serves as a relative measure of 
the degree of graphitization, that is, the degree to which the filler crystal structure 
approaches that of a theoretically perfect crystal (dc for a perfectly graphitized crystal 
is 3.354 A). The interlayer spacing, measured by X-ray diffraction, was taken from 
reference 3. 


Thermal Conductivity 

Accurate measurement of the thermal conductivity k of graphite is both difficult 
and tedious. Reference 9 points out, however, that the thermal conductivity of a graphite 
is uniquely related to its electrical conductivity Og, and therefore can be predicted from 
it. This is a particularly valuable situation because accurate electrical conductivity data 
are relatively easy to obtain. Accordingly, the thermal conductivities of the present 
graphites were determined from experimental measurements of electrical conductivity. 
The equation used to relate these two parameters is given by 

k = 1.513 X 10-3 _ 0.1388 X 10’6 

and was specifically developed for the present graphites. Its development is given in 
appendix A. All measurements of electrical conductivity were made in the across -grain 
direction by a four -probe direct -current technique and were performed by the Oak Ridge 
National Laboratory on solid cylindrical samples furnished by NASA. This four -probe 
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technique uses direct -current leads brought into contact with the end faces of the cylindri- 
cal sample and knife-edges of known separation contacting the surface of the sample from 
which the potential drop along the specimen axis is obtained. Electrical conductivity mea- 
surements were made on three samples of each graphite and the results averaged. 

ABLATION -PERFORMANCE MEASUREMENTS 
Test Environment 

The test environment consisted of a nominal Mach 2 airstream, with a nominal total 
enthalpy of 2.17 MJ/kg and a stagnation pressure of 5.6 atm (5.67 x lO^ N/m2). Heating 
rates were nominally 513 W/cm2. This condition was obtained in the arc -heated materials 
jet at the Langley Research Center which is described in reference 10. Stagnation pres- 
sures and heating rates in the test stream were measured with a calibrated pressure 
probe and a calibrated Garden foil asymptotic calorimeter, respectively. Stream enthal- 
pies were calculated from the pressure and heating rate according to the theory of Fay 
and Riddell (ref. 11). 


Test Specimens 

The artificial graphite test specimens were machined from graphite billets as sup- 
plied by the manufacturer. Each specimen consisted of a 1.270 -cm -diameter hemisphere- 
cylinder, 1.905 cm long, with a 0.635-cm nose radius. During an ablation test, the graph- 
ite test specimen was mounted in a 1.905 -cm -long phenolic -asbestos insulator which was, 
in turn, mounted in a water-cooled holder (see fig. 1). All graphites were mounted so 
that the test stream impinged on the specimen normal to the preferred orientation of the 
basal plane surfaces (i.e., in the across-grain direction), although reference 2 demon- 
strates that this attention to orientation is not really necessary provided, as in the present 
case, that there is no efficient path for thermal conduction from the specimens. Only one 
ablation-performance test was made on each graphite. 

Instrumentation and Data Analysis 

The instrumentation used for measuring ablation response of the graphite specimens 
consisted of a bench micrometer (direct reading in inches to 0.0001 inch (0.00025 cm)), a 
motion-picture camera with a framing rate of 200 pictures per second, and a continuous 
recording photographic pyrometer, the theory and principle of which are described in ref- 
erence 12. The parameters selected for evaluating ablation performance were total 
length change A.I resulting from a 30 -second exposure to the test environment, linear 
stagnation-point recession rate s, steady-state stagnation -point surface temperature T, 
surface texture t, and degree of gouging F . 
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Total length change was obtained from length measurements of the specimen before 
and after test using the bench micrometer. Linear recession rate was determined from 
the motion-picture -film record of the eroding specimen as follows. Specimen length as 
a function of time was obtained from the film records of the specimen with the aid of a 
motion analyzer. These data were then plotted as shown in figure 2. Customarily, spec- 
imen length initially increased because of thermal expansion of the graphite, but soon 
passed through a maximum and thereafter decreased continuously because of erosion. 
Recession rate was taken to be the slope of the curve over the linear portion of the length 
decrease with time as determined by a least-squares fit of the data. 

Steady -state stagnation -point surface temperature was determined from the motion- 
picture -film record taken with the photographic pyrometer. Stagnation -point surface 
temperatures as a function of time were obtained from this record and plotted as shown 
in figure 2. For determining these temperatures, the emissivity of graphite was taken 
to be 0.94 (refs. 13 and 14). The representative stagnation -point surface temperature 
was taken as the steady -state maximum temperature attained during the run. In all cases, 
this temperature was reached during the linear portion of the length -time curve. 

Surface texture is a measure of the degree of small-scale unevenness or roughness 
of the surfaces of the after -test specimens, as assigned by visual inspection. Numbers 
range from 1 (smooth) to 3 (rough). Degree of gouging is a measure of the extent to which 
the after-test specimens showed gross material loss from their surfaces, as evidenced 
by gouges, grooves, or missing chunks of material. Numbers from 1 (no gouging) to 3 
(severe gouging) were assigned by visual inspection. The numerical designations for 
both surface texture and degree of gouging were assigned in accord with the photograph 
of the representative after-test specimens shown in figure 3. 

Tsqjical Test Sequence 

A typical experimental test sequence proceeded as follows. The magnetically 
rotated electric-arc air heater was started and 12 seconds were allowed for the arc to 
stabilize and the test-stream environment to become established. The motion -picture 
camera and photographic pyrometer were started and the test -stream heating rate was 
measured. The graphite specimen to be tested was inserted into the stream and allowed 
to remain for 30 seconds. After retraction of the specimen, a final heating-rate mea- 
surement was made. Each event in this sequence was executed automatically and con- 
trolled by a preset programmer. 

EXPERIMENTAL RESULTS AND DATA ANALYSIS 

The results of the material -property measurements and ablation-performance mea- 
surements are shown in table I. All data are complete on the first 40 graphites; however, 
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the last five were added during the experimental program and complete material -property 
data have not been obtained on them. Accordingly, in the following discussions regarding 
correlations of material properties and ablation performance, only the data on the first 
40 graphites are considered. 

In addition to the ablation test data shown in table I for the two Graph -i-tite "G" 
grades, separate tests were made to compare their ablation performance in more detail. 
These additional ablation-performance data are shown in table II and discussed in 
appendix B. 


Material-Property Correlations 

It is apparent from the description of the material -property measurements in the 
preceding section that, although conventional definitions of the properties were used, 
many of the properties ought to be closely related (i.e., such properties as pore area Ap 
and surface area S, the various alternate definitions of a characteristic pore dimension, 
etc.) Accordingly, it was deemed an important part of the present study to determine, 
where possible, correlations of these material properties, thus reducing the total number 
of property variables requiring consideration. In order to provide the maximum useful- 
ness from this reduction, correlations were sought predicting the more -difficult -to - 
measure properties in terms of the simpler -to -measure properties. For instance, since 
the surface area S is simpler to measure than the pore area Ap (because only the 
first part of the adsorption isotherm is required to determine S), it is desirable to 
develop a correlation predicting Ap in terms of S. For simplicity, correlations were 
sought predictii^ each property in terms of only one other prcperty instead of a combina- 
tion of other properties. Also, since plots of each and every pair of properties demon- 
strated that, in all cases, where a relationship existed it was approximately linear, only 
linear correlations were considered. In developing such correlations, the method of least 
squares was used. The correlations developed are given by the following equations: 


Prj, = 4.489 - 0.664dc 

R= 0.01 

(1) 


^-0.1066 + 0.00924e 

R = 10.6 

(2) 


-0.00142 + 0.006540^.^ 

R= 7.6 

(3) 


0.7970 - 0.3980pp 

R = 11.1 

(4) 

Vp = -0.0046 + 0.0029 5S 

R = 12.5 

(5) 

) 

"0.072 + 0.846S 

R= 4.7 

(6) 

Ap - ^ 



^.37 + 238.0vp 

R = 10.8 

(7) 
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- 0.2 + 2.21rj^ 

R= 4.2 

(8) 

|l2.10 + O.67O0^ij. 

1 5^1 

11 

CO 

• 

00 

(9) 

< 10.41 + O.71O0Hg 

R = 5.7 

(10) 

(98.29 - 43.36pp 

11 

o 

« 

00 

(11) 

, = 120.68 - eo.oipp 

R= 7.2 

(12) 


pLlO.43 - 53.15pp 

R = 10.2 

(13) 


3.49 + 0.889e^jj. 

R= 6.6 

(14) 

1 

! 

^.71 + 135.7Vp 

R= 2.2 

(15) 


Predicted values from these correlations are compared with measured values in fig- 
ure 4. All correlations have a correlation coefficient greater than 0.92 except the cor- 
relation of pore area in terms of pore volume which has a correlation coefficient of 0.89. 

It is pointed out that a correlation was developed relating to d^ (eq. (1)) even 

though prj, is defined in terms of d^. because such a linear expression might be desir- 
able for some purposes; the remaining correlations were developed for pairs of variables 
having no previously known functional relationship. It is noted that four material proper- 
ties (Vp, Ap, e, and 0jjg^ can be correlated satisfactorily by two or more expressions. 
It is also interesting to note that of the eight material properties correlated, four (Vp, e, 
®air> ®Hg) correlated satisfactorily in terms of the bulk density pp. It is 

not surprising, then, that bulk density should have become a most important property in 
selecting graphites for consideration as ablators. 

Since two variables are obviously not independent if one can be predicted from the 
other, the total number of properties under consideration for describir^ the various 
artificial graphites can be reduced by using the correlations shown in equations (1) to (15). 
Hence it is seen that out of the original 18 material properties only the 10 properties 
pp, Rp, Gm, S, r^n, rp. A, d^, and k are truly independent. This reduc- 

tion in material properties requiring consideration represents a savings in both experi- 
mental and computational time. 

Treatment of Ablation -Performance Data 

In order to determine which of these 10 material properties have the most pro- 
nounced effect on graphite ablation performance in the present test environment, attempts 
were made to relate the five ablation-performance parameters AZ, T, s, t, and T 
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(considered as dependent variables) to the 10 material properties (considered as inde- 
pendent variables). Standard methods of stepwise multiple, curvilinear regression were 
used (ref. 15). Each material property was considered in the regression equation to the 

first power alone, and to both the first and second power according to the general exp res - 
10 

Sion y = c + ^ (^i^i ^i^i^)> where y represents an ablation performance parameter, 

i=l 

xj the ith material property being considered, and aj, bj, and c the coefficients 
calculated from the regression. Although this form of equation obviously cannot yield 
information regarding a true functional relationship or mechanism, it can identify those 
variables which are functionally related in a statistical sense, and, furthermore, can 
indicate the direction of the functional dependence - that is, whether an increase (or 
decrease) in the value of a certain material property will result in an increase or decrease 
in the value of an ablation parameter of interest. The main purpose in formulating the 
regression equations was not to determine a true functional mechanism nor to predict 
ablation performance (though it can, in some cases, be useful for this), but was to identify 
those material -property variables which have the strongest effect on ablation performance. 

Since the computations required are extensive, all computations were performed on 
a digital computer using a modified form of a computer program written by Efroymson 
(ref. 16). One very useful feature of this computer program is that it automatically tests 
against a preselected significance level the significance for entry of each of the x^ and 
Xj^^ into the regression equation. This insures that all of the Xj and x^^ the fi nal 
regression equation are statistically significant at the preselected level. Accordingly, 
all variables having only a random effect (at the selected significance level) are auto- 
matically excluded. The computer program used is given in appendix C. 

RESULTS AND DISCUSSION 

Dependence of Ablation Performance on Material Properties 

Regression equations for the five ablation parameters were formulated in terms of 
the 10 independent property variables at a significance level of 1 percent. A 1 -percent 
significance level means that if a given property appears in the final regression equation, 
there is only a 1 -percent statistical probability of its having been erroneously included, 
or, in other words, there is only one chance in 100 of that variable having but a random 
effect on the dependent variable (i.e., on the ablation parameter). It is important to 
recognize, however, that such a significance criterion does not guarantee that every vari- 
able which has a true effect on ablation performance has been discovered; neither does it 
absolutely guarantee that every variable which is singled out by the significance test as 
having an effect does, in fact, have such an effect. What it does guarantee is that, for 
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each of those variables singled out, there is only one chance in 100 (1 -percent statistical 
probability) of that variable not having a real effect. 

Obviously a significance level other than the present 1 -percent level could have 
been chosen; the precise choice must be made on the basis of the degree of risk one is 
willing to take in incorrectly assuming that a variable has an effect when it really does 
not. The reasoning behind the present choice of 1 percent was that if a graphite manu- 
facturer were to attempt to fabricate an improved graphite by accepting the property 
variables singled out by the regression equations as important, and were to proceed by 
adjusting his formulation and fabrication techniques to optimize these variables, he might 
reasonably desire the assurance that there be only one chance in 100 that he is proceed- 
ing in an unprofitable direction. A matter of some interest, however, is that even though 
the present regression equations were formulated at a significance level of 1 percent, it 
turns out that all variables are also significant at the 0.5 -percent level. In other words, 
if one accepts the singled-out property variables as being important, there is actually 
only one change in 200 of his being wrong. 

Regression equations were formulated for all five ablation parameters considering 
the material -property variables to the first power only, and also to both the first and 
second powers. In both cases, the material properties singled out as significant were 
identical; therefore, only those regression equations obtained by considering the material - 
property variables to the first power alone are presented. The resulting regression 
equations, with their average absolute percent residuals R, are given by 


Al = 2.073 - 0.98 5pp + 0.437Gm 

R = 12.3 

(16) 

T = 2169.9 + 157.8Gm + 0.0201A - 67.4k 

R = 1.5 

(17) 

s = 0.05156 - 0.02145pp + 0.01596Gm 

R = 10.7 

(18) 

T = 0.7 + 2.7Gxxi 

R = 14.9 

(19) 

r = 6.2 - 3.1pp + 2.5Gm 

R = 28.0 

(20) 


Equations (16) to (20) clearly reveal the importance of maximum grain size Gm 
and bulk density pp in determinir^ ablation performance. Maximum grain size is the 
only property appearii^ in all five regression equations. As Gm increases, so do AZ, 
T, s, T, and T. This result is in agreement with the currently accepted practice of 
selecting fine-grained graphites for ablation materials. Equation (19) indicates that sur- 
face texture t is a function of Gm only, which certainly appears reasonable. That 
Ai and s are strongly dependent on Gm (see eqs. (16) and (18)) could imply, among 
other things, a contribution of particle removal to the ablation process (ref. 17). Bulk 
density pp also affects AZ and s, with equations (16) and (18) indicating that an 
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increase in pp produces a correspondir^ decrease in Al and s. This is most rea- 
sonable because for a given rate of reaction (rate of removal of mass) dimensional 
changes will be smaller the denser the graphite. This predicted behavior of higher den- 
sities resulting in smaller values of Al and s is consistent with current practice of 
selecting high-density graphites for ablation applications. Also, as seen from equa- 
tion (20), higher densities tend to reduce the degree of gougii^ of an ablatii^ graphite. 

Ash content A is seen from equation (17) to affect surface temperature T, with 
higher values of A resulting in higher surface temperatures. This result is to be 
expected, since the presence of chemical impurities in an artificial graphite is known to 
catalytically accelerate its oxidation rate (refs. 1 and 2), driving the surface temperature 
higher by the exothermicity of the reaction. A reasonable upper limit on ash content can 
be derived from equation (17). For instance, since an ash content of 100 ppm produces 
an increase in surface temperature of only about 2 K, it seems reasonable that one might 
want to maintain the ash content at or below that level, which is not at all difficult to 
achieve. It is also seen from equation (17) that thermal conductivity k has a negative 
effect on surface temperature. This result, again, is most reasonable since higher 
values of k will serve to increase the rate of heat transfer from the stagnation point 
to the sides and other parts of the graphite specimen. Thus, on the whole, it appears 
that the material properties primarily responsible for determinii^ the ablation perform- 
ance of a t5q)ical artificial graphite are Gm, Pp, A, and k. 

Since artificial graphites which develop roughened surfaces or become gouged dur- 
ing ablation are often undesirable in practical applications, criteria by which such graph- 
ites can be eliminated from consideration are valuable. Equations (19) and (20) provide 
such criteria. For instance, if graphite surface texture is required to be t ^ 1.3^ appli- 
cation of equation (19) shows that all those graphites with grain sizes Gm > 0.22 mm 
should be eliminated. Examination of table I shows that there are 11 graphites of the 
first 40 having Gm > 0.22 mm. Of these 11, only two had smooth surfaces after testing; 
at the same time, only one graphite (4007) would not be eliminated which had a slightly 
roughened surface (t = 2). Now, in addition, if it is also required that the degree of 
gouging be r s 1.5, application of equation (20) indicates that additional graphites should 
be eliminated. Since F is a function of both density Pp and maximum grain size Gm, 
the criteria for elimination of graphites will be in terms of both variables. The simplest 
way to establish these criteria is to calculate, for each different value of maximum grain 
size, that value of density below which if the actual graphite density falls, the graphite is 
to be eliminated. A list of such density values for the 11 different values of maximum 
grain size of the present graphites is shown in table HI. In addition to the 11 graphites 
already eliminated by the requirement that t S 1.3 (Gm = 0.22 mm), eight more will be 
eliminated by the added requirement that F ^ 1.5 (given by the criteria in table m). The 
one graphite (4007) not eliminated previously which had a slightly roughened surface is 
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now eliminated by this additional requirement. Note also that this graphite is severely 
gouged (r = 3). Considered by itself, the requirement that r ^ 1.5 eliminates all but 
one graphite which should be eliminated (H205-85 with r = 2), and eliminates only three 
graphites which were, in fact, not gouged. 

Elimination of the 19 graphites indicated by these two requirements on surface 
texture and degree of gouging produces a subset of the original data set which is com- 
prised of graphites virtually all of which erode with smooth, ungouged surfaces. This 
subset may then be analyzed in a maimer similar to that of the original data set — that 
is, regression equations can be formulated for the ablation parameters AZ, T, and s, 
in terms of the 10 material properties. This was done, as before, at the 1 -percent sig- 
nificance level, with the results that Al is solely a function of Pp, that T is a func- 
tion of rm and A, and that s does not correlate with any material property. How- 
ever, because the reduced number of data sets for this analysis (21 as opposed to the 
original 40) makes it more likely that the high significance level (1 percent) will cause 
important variables to be omitted from the regression equations, this significance level 
was relaxed somewhat to the 2. 5 -percent level (the next lower significance level at which 
statistical tests are customarily made). At this 2.5-percent level, the following regres- 
sion equations were obtained: 


Al = 1.587 - 0.583pp + 1.82 x IQ-^A - 2.97 X 10'^ rm 

R = 4.7 

(21) 

T = 2843.1 - 5.475rm + 0.0294A - 253.4pp 

R = 1.0 

(22) 

s = 0.03457 - 0.01028pp 

R = 5.5 

(23) 


Bulk density Pp is seen to affect AZ and s as before; that is, the greater the 
density, the smaller AZ and s. Density is also seen to influence T, with lower values 
of T the higher the density. Also as before, an increase in ash content A is seen to 
produce a higher surface temperature; but A is also indicated as having an effect on Af 
(higher A producing larger Al). This is, of course, again consistent with the known 
effects of chemical impurities catalytically accelerating oxidation rate (refs. 1 and 2). 

The reason why A is not also indicated as having an effect on s is very likely asso- 
ciated with the fact that s, being a much more difficult parameter to measure than Al, 
is known less accurately, and as the error in a dependent variable increases, the sensi- 
tivity of the regression analysis decreases. In fact, it so happens that if the significance 
level criterion is relaxed to about the 7.5-percent level, A is the next property variable 
to enter the regression equation for s, and, as should be the case, an increase in s is 
predicted by an increase in A. The influence of rm on Al and T is interesting and 
bears further investigation. However, it is gratifyingly assuring that a change in rm 
causes both Al and T to change in the same direction, since a change in opposite 
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directions for AZ and T is most unlikely. Hence, taken together, equations (21) 
to (23) reveal that Pp, A, and rm are the material properties of importance in deter- 
mining the ablation performance of a t 5 rpical artificial graphite intentionally selected to 
ablate with a smooth, ungouged surface. 

In addition to the previously discussed theoretical limitations concerning regression 
equations in general, there are other possible, more practical limitations on the present 
regression equations (eqs. (16) to (23)). Primary among these is the possibility that not 
all of the important material properties may have been included in the list of properties 
considered. Several potentially important properties not included are: the raw materials, 
filler -to -binder ratio, grain-size distribution (instead of merely the maximum grain size), 
the distribution of impurities between filler and binder carbons, the lattice dimensions of 
the binder, the nature of the bonding between the filler and binder, the type of porosity, 
the gas permeability, tensile strength, shear strength. Young's modulus, etc. Some of 
this information is difficult or, in the case of proprietary information, impossible to 
obtain. Even if all this information were available, some is difficult to specify quantita- 
tively in meaningful terms (for instance, the distribution of impurities between filler and 
binder carbons). But, in spite of the possible omission of some potentially important 
material properties, the present analysis has revealed five material properties to be 
important in determining ablation performance (Ghij Pp> A, k, and rj^^. Of par- 
ticular importance, also, is that the present method of analysis has been shown to be both 
powerful and efficient, at no time contradicting either experience or reason in predicting 
a given material property to be important or in predicting the direction of the effect 
which a given material property has on ablation performance. 

The information gained from the regression equations in this study should be useful 
in selecting artificial graphites for consideration for use in environments similar to that 
of the present test environment, and also as a guide to future research. Certainly the 
present method of analysis appears particularly effective. For a more complete picture 
of the various material properties affecting graphite ablation performance, it would be 
desirable to apply the present approach to ablation data obtained in different environments 
spanning a wide range of pressures and enthalpies. Similar regression equations could 
be developed for these different environments and, as in the present study, criteria could 
be established for the elimination of those graphites which are likely to develop unde- 
sirably rough surface, unacceptably high degrees of gouging, or other undesirable charac- 
teristics. A more fundamental investigation could perhaps be conducted by carefully con- 
trolling the fabrication of different graphites, with each graphite differing systematically 
in its material properties. The severe difficulty to be encountered in this approach, 
however, is that it is almost impossible to alter the value of any one material property 
without at the same time altering many others as well. A combination of these two 
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approaches - controlled fabrication and statistical treatment of the data - may well be 
the most promising approach for future work. 


Comparisons With Previous Work 

Several other investigations also attempting to determine the effect of material 
properties on the erosion behavior of graphite in high temperature dynamic environments 
have been reported (refs. 18 to 20). Slosarik and Swope (ref. 18) and Chase (ref. 19) 
investigated the effects of certain material properties on rocket -nozzle erosion and 
McVey, Auerbach, and McBride (ref. 20) investigated material -property effects on abla- 
tion. Slosarik and Swope considered the properties of density, porosity, pore-size dis- 
tribution, tensile strength, and impurity content, but could find no direct, unambiguous 
effects. However, they concluded from the trends of their data that high density, low 
porosity, low impurity content, high tensile strength, and a small number of larger sized 
pores are important for good erosion resistance. Chase (ref. 19) considered flexural 
strength, density, and porosity and concluded that all three "play their part in determin- 
ing erosion resistance." McVey, Auerbach, and McBride (ref. 20) considered the effect 
of graphite microstructure on ablation performance, and concluded that grain size, filler - 
to-binder ratio, porosity, porosity type, grain orientation, and degree of graphitization 
all affect ablation performance. 

All three of these reported investigations suffer from the weakness that in attempt- 
ing to determine which material properties most strongly affect ablation performance, 
each material property was considered separately, one by one, for its individual effect 
on ablation performance. In many cases, the data presented are ambiguous and open to 
several different interpretations. Consequently, many of the conclusions reached, though 
sufficiently reasonable, are actually based on theoretical or intuitive arguments and lack 
the necessary supporting experimental data. Two of the conclusions reached by McVey, 
Auerbach, and McBride (ref. 20), however, conflict with results obtained in the present 
study or with those of other studies. These conclusions are that grain orientation is 
important, and that the greater the degree of graphitization, the more readily the graphite 
will ablate. In conflict with their first conclusion are the experimental results of refer- 
ence 2 demonstrating that grain orientation has little effect on the ablation performance 
of artificial graphite. Conflicting with their second conclusion is the discussion in ref- 
erence 3 which points out that the oxidation rate of a graphite should, if an3rthing, decrease 
with increasing degree of graphitization. Furthermore, since degree of graphitization is 
directly related to interlayer spacing dc (ref. 3) which does not appear in the regres- 
sion equations of the present study, it would appear that degree of graphitization is not 
one of the more important material properties affecting ablation performance. One pos- 
sible explanation why McVey, Auerbach, and McBride arrived at these two debatable con- 
clusions is that they studied an extremely broad diversity of carbonaceous materials - 
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commercial graphite, a carbon felt impregnated with pyrolsdic carbon, and a composite 
material consisting of carbon fibers in a graphitized matrix - as opposed to studyii^ a 
single class of carbonaceous materials such as commercial artificial graphite, for 
instance. For the materials they studied, the precise significance of a number of the 
material properties is vague. Certainly the conventional concept of grain orientation, 
and perhaps even degree of graphitization lose much of their significance.. 

It is of some interest to note that all three investigations cited above (refs. 18 
to 20) report porosity to be an important variable, with the first two investigations also 
reporting density to be important. This is not surprising considerii^ that it has been 
shown in the present study that total porosity e can be predicted in terms of bulk den- 
sity Pp to less than 1 percent and open porosity (^air ^Hg) predicted to 

within 10 percent or less (eqs. (11) to (13)). Hence, the conclusion that porosity is impor- 
tant is merely an alternative way of stating that density is important. The two conclu- 
sions should therefore not be considered as separate — at least, not until a separation 
can be made on the fundamental basis of mechanism. 

Predictions of Ablation Performance of Several Graphites 

As mentioned previously, several graphites were added during the experimental 
program and complete material -property data were not determined on them. These 
graphites, the last five in table I, were included in the investigation because of special 
characteristics of possible value: AGOT graphite was chosen because it is reported to 
be unusually resistant to the propagation of cracks caused by thermal shocking (ref. 21); 
AXF-5Q graphite was chosen because of its exceedingly fine grain size; and the three 
Graph-i-tite graphites were chosen because of reports (unpublished) of high erosion 
resistance in ablation tests. Also as mentioned previously, since the property data on 
these graphites are incomplete, they were not included in the development of the regres- 
sion equations (eqs. (16) to (23)). Nonetheless, some property data are available on these 
graphites (table I) , and even though the expressed purpose in developing the present 
regression equations was simply to isolate those material properties having the strongest 
effect on ablation performance, it is possible to predict ablation performance for some 
of these graphites and to compare this predicted performance with actual measurements. 

Ablation-performance predictions can be made for AGOT, AXF-5Q, Graph-i-tite ”A,” 
and the Graph-i-tite "G" graphite which has a maximum grain size of 0.840 mm. (The 
Graph-i-tite "G" graphite with a maximum grain size of 0.203 mm was tested in a differ- 
ent series of tests at a somewhat higher heatii^ rate than the other graphites (562 W/cm2) 
and therefore predictions of ablation performance are not strictly comparable.) Ablation 
performance predictions have been made using equations (16) to (20), where appropriate, 
with the results shown in table IV. The only acceptable agreement between predicted and 
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measured ablation performance is for surface texture t and degree of gougii^ r. 
Agreement for surface temperature T is not bad; but, then, since the spread in test 
temperatures is not large, the predicted temperatures cannot be too far off. For AZ 
and s, the predicted values for the large -grained graphites are too large and the pre- 
dicted value for the small-grained graphite is too small. This lack of agreement is pos- 
sibly related to the special nature of these graphites, with the direction of the discrepancy 
undoubtedly being due to the strong dependence of AZ and s on Gm- 

The application of the criteria established by the requirements on surface texture 
and degree of gouging discussed previously (t S 1.3 and F S 1.5) would eliminate AGOT, 
Graph-i-tite "A," and the large grained form of Graph-i-tite "G" from consideration. 

This is as it should be, since these graphites developed roughened, gouged surfaces during 
ablation. The remaining two graphites (AXF-5Q and the fine grained form of Graph-i-tite 
"G") both have smooth, ungouged surfaces. Equations (21) to (23), then, ought to be appro- 
priate for predicting their ablation performance. This can be done only for AXF-5Q, 
however, because, as previously stated, these equations are not strictly applicable to this 
fine-grained form of Graph-i-tite ”G" since it was ablation tested at a heating rate differ- 
ent from that at which the regression equations were formulated. In attempting these 
predictions, it is noted that mean pore radius r^ appears in the equations for both AZ 
and T, and since r^^ is not known for AXF-5Q, AZ and T cannot be calculated. 
However, s is a function of pp alone (eq. (23)) and, therefore, values of s can be 
calculated. Its predicted value, given in table IV, is seen to be lower than the experi- 
mental value, but not as low as that predicted by equation (18), which includes a term for 
the maximum grain size. It is worth reemphasizing, however, that the main value of the 
present regression equations is in identifying the important material properties and the 
directions in which they affect ablation performance, and not in predictii^ ablation per- 
formance itself. 


CONCLUSIONS 

Measurements of 18 material properties and of ablation performance in one test 
environment have been made on 45 commercial, artificial graphites. As a result of these 
measurements and an analysis of the data, the followii^ conclusions and recommenda- 
tions are drawn. 

1. Correlations between pairs of the 18 graphite material properties have been 
developed where possible, thus enabling certain properties to be predicted in terms of 
others. The original 18 material properties were thereby reduced to 10 mutually inde- 
pendent properties. 

2. A correlation based on 10 selected graphites was developed relating thermal 
conductivity to electrical conductivity. This correlation is in good agreement with 
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existing correlations in the literature, but was formulated specifically for the present 
graphites which have thermal and electrical conductivities extending to lower values than 
had been included in previous correlations of thermal and electrical conductivity. 

3. Regression equations have been developed at the 1 -percent significance level 
relating the ablation performance of artificial graphite in the present test environment 
(Mach 2 airstream, stagnation pressure of 5.6 atm, and enthalpy of 2.2 Mj/1^) to graphite 
material properties. These regression equations which, in fact, are even significant at 
the 0.5-percent level, demonstrate that ablation performance depends on maximum grain 
size, bulk density, ash content, and thermal conductivity. 

4. Criteria have been established from these regression equations for eliminating 
those graphites likely to develop rough surfaces or become gouged during ablation. After 
elimination of these graphites, a second set of regression equations have been developed 
relating ablation performance to material properties. These equations reveal that abla- 
tion performance of artificial graphites intentionally selected to ablate with smooth, 
ungouged surfaces depends on bulk density, ash content, and mean pore radius. 

5. Maximum grain size is an important property in determining ablation perfor- 
mance, with an increase in maximum grain size causing an increase in total length change, 
surface temperature, recession rate, surface texture, and degree of gouging. To ensure 

a smooth surface, maximum grain size should be kept below 0.22 mm. 

6. Bulk density is also of prime importance in determining ablation performance, 
with higher densities resulting in lower values of total length change, recession rate, and 
degree of gouging. 

7. High ash content produces higher surface temperatures, while high values of 
thermal conductivity and mean pore radius tend to produce lower surface temperatures. 

8. In spite of certain theoretical and practical limitations, the present method of 
analysis (stepwise multiple regression) for identifying the more important graphite mate- 
rial properties affecting ablation performance is both powerful and efficient. Its sound- 
ness is demonstrated by the fact that it properly predicts ablation performance to be 
related to certain material properties known from theoretical or experimental considera- 
tions to affect ablation performance. 

9. In order to determine other material properties which might also have an impor- 
tant effect on ablation performance, a more extensive list of properties than the 18 con- 
sidered here could be employed. Also, testing should be conducted in a wider range of 
environments. A desirable addition to such a study would be the systematic, controlled 
fabrication of different graphites, differing, insofar as possible, in only one property at 

a time. 
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10. Graph-i-tite "G" graphite is particularly erosion resistant, even more so 
than the usual standard of comparison, ATJ graphite. Whereas the larger -grained 
Graph-i-tite "G" (maximum grain size of 0.840 mm) became rough during ablation, the 
finer-grained Graph-i-tite "G" (maximum grain size of 0.203 mm) remained smooth, 
thus making it particularly attractive for future study. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 7, 1972. 
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APPENDIX A 


PREDICTIONS OF THERMAL CONDUCTIVITY FROM 
ELECTRICAL CONDUCTIVITY 

o 

Reference 9 reports the relation k = 1.297 X 10 Oj, for predicting thermal con- 
ductivity from electrical conductivity (at room temperature). In order to verify this 
relation for the present graphites, both thermal conductivity and electrical conductivity 
were experimentally measured on 10 of the graphites. These 10 graphites were selected 
from the total number of graphites in the present study to span the full range of measured 
electrical conductivities in approximately equal increments. Thermal-conductivity mea- 
surements were made using a comparative longitudinal heat -flow apparatus and were 
performed by the Oak Ridge National Laboratory on the same samples used for the 
electrical-conductivity measurements. This apparatus consists of a vacuum chamber 
containing an axial column of two instrumented Armco iron heat meter bars between 
which the specimen is compressed. A heater on one of the meter bars induces heat flow 
in the column and axial temperature distributions are measured with thermocouples on 
the bars. 

The resulting thermal and electrical-conductivity data are plotted in figure 5, with 
the curve through the data beir^ given by the expression 

k = 1.513 X 10-3 _ 0.1388 x 10"® 

This equation was developed from the data by the method of curvilinear least squares. 

It was formulated specifically for the present graphites because they include thermal and 
electrical conductivities extending to lower values than had been included in previous cor- 
relations of thermal and electrical conductivity. Predictions of thermal conductivity 
from this equation compare favorably with predictions from the equation of reference 9, 
and even more favorably with predictions from an equation developed by the Oak Ri<^e 
National Laboratory 

k = 1.56 X 10"3 (Tj, - 0.266 X 10~6 

This latter equation was developed from more extensive data, including data on the 10 
NASA graphites, as well as on other graphites with electrical conductivities extending as 
high as 1940 (S2-cm)"^. 
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APPENDIX B 


ABLATION PERFORMANCE OF Graph-i-tite "G" 

During the course of the present investigation it was observed that the Graph-i-tite 
"G" graphite being tested (maximum grain size 0.840 mm) was particularly more erosion 
resistant than the other graphites under test. However, anomalously, it also eroded to 
produce a very rough surface. It was reasoned that its rough surface must be due to its 
large grain size, but that its high erosion resistance must be due to other factors. Since 
Graph-i-tite "G" is also available with the smaller maximum grain size of 0.203 mm, it 
was further reasoned that this fine-grained form might be as erosion resistant as the 
large-grained form, but at the same time erode with a smooth surface. Accordingly, 
some of the fine-grained material was obtained, and a brief test series conducted in 
which the ablation performances of Graph-i-tite "G" (Gm = 0.840 mm), Graph-i-tite "G" 
(Gm = 0.203 mm), and ATJ were compared. (The ATJ was selected because it is a com- 
monly used standard of comparison.) Nominal heating rate (562 W/cm2) was somewhat 
higher in these tests than in the others. 

Results of this test are shown in table H, where it is observed that, within experi- 
mental scatter of the test, the recession rates s of the two Graph-i-tite "G" graphites 
are virtually identical, but significantly lower than that of ATJ. Also, the total reces- 
sions A2 of the two Graph-i-tite *'G" materials are close, but, again, are much lower 
than that of ATJ. It appears from these tests, therefore, that the fine-grained form of 
Graph-i-tite "G,” which, in fact, did erode with a smooth, ui^ouged surface (see table I), 
may offer improved ablation performance compared with that of other graphites cur- 
rently in use. 
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APPENDIX C 


COMPUTER PROGRAM FOR STEPWISE MULTIPLE 
CURVILINEAR REGRESSION 

The computer program used is a modified form of that written by Efroymson 
(ref. 16). The computational procedure is the same as that originally written by 
Efroymson, but a number of options have been added to the program making it more 
flexible. 

Certain important details of the program are that the number of independent vari- 
ables is limited to 49, and the number of sets of observations is limited to 100. A weight- 
ing factor can be assigned to each set of data if desired. A significance level must be 
chosen and a corresponding F -level (see, for instance, ref. 15) must be set for the entry 
of a variable into the regression equation. The computation then considers each inde- 
pendent variable step by step with a variable entering the equation only if it produces a 
statistically significant reduction in the residual sum of squares about regression (at the 
chosen F -level). 

The data are read in according to the format 6F12.5, with all the independent vari- 
ables of a set read in first, followed by the dependent variable, and followed, in turn, by 
the weighting factor, if any. The control card is written as follows: 


Col. no. 

Variable 

Description 

1-10 

TOL 

Tolerance (set at 0.001) 

11-20 

EFIN 

F -level to enter variables into regression 

21-30 

EFOUT 

F -level to remove variables 

31-35 

NOPROB 

Problem number 

36-40 

INVAR 

1 greater than the number of independent variables read in as 
data 

41-45 

NODATA 

Number of data sets 

48 

IFWT 

1, if no weighting factors assigned 

50 

IFSTEP 

1, if intermediate steps not to be printed 

52 

IFRAW 

1, if residual sums and squares not to be printed 

54 

IFAVE 

1, if average not to be printed 

56 

IFRESD 

1, if residual squares and cross products hot to be printed 

58 

IFCOEN 

1, if partial correlations not to be printed 

60 

IFPRED 

1, if predicted values not to be printed 

62 

IFCNST 

1, if constant term in regression equation is assumed to be zero 

64 

IFDATA 

1, if input data not to be printed 
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APPENDIX C - Continued 

Variable De s cription 

IFPWR 1, if second power of variables not to be considered 

IFPLOT 1, if predicted values not to be plotted against measured values 

The F-level is based on 1 and (NODATA - NOVAR) degrees of freedom where 
NODATA is the number of data sets and NOVAR is 1 greater than the total number of 
independent variables considered in the regression equation. (In the case where the sec- 
ond power of the variables is to be considered, the number of independent variables con- 
sidered in the regression equation is obviously twice the number of variables appearing 
as raw data.) The F-level for entering variables must be greater than or equal to the 
F-level for removing variables. A complete listing of the computer program is given 
below. 

C SIMPLF STFP'VISF mult rPLE-^UPVlL I NFAR RFGRFFSION 

r CONSIDFRf FACH VARIApL^^ to 1 ft P^'r'FR, ^NO IF DESIRFD 

D r MFNS ion data (EO ) , Vp*CTOR CFO ,50 ) ^ AVF ( FO ) ^ SI GMA (Fn ) , 

1 COFN (FO),SIGMrO(Fn),TNDf^X(F^)iDATAS(100,FO),X(100)^Y(10n5 
DIMFMSIOM TN(?> 

r INVAP = 1 + NO. OF I Mr>,c-pprKjnpi^ y VARIAPL^^S RfAD in A^ data 

r NOVAR = 1 + TOTAL NO. Qf IMOFPFNDFNT VARIARLi^^S CONS I DFRFD by PROGRAM 

r MOOATA = NF. OF OATA fptS 

C FFIN AND f^f^OUT ARp^ f/ssf^j qN 1 AND (NODATA - NOVAR) DFGR'^f^S OF FREEDOM 

inn RFAD TOL 1 Ef tn , r-FOUT , NOf^^ROb , I NVAR , MOO at a , 

1 IFWT, IFST^n, IPRAW. IEAVE» I FRF SR, IFCOENf IFPRFD, IFCNST, IFRATA. iFPWRf 

F j PPLOT 

IF (FOf.F) QQD.Oon 
’^9P ST"^P 
PRO CONT^NU^ 


r: 

IFM/T= 1 4 

T 

HFN 

HIHTF 

= 1 . 










r 

I fst^p= 

1 

♦ no 

NOT 

PR I M 

T 

^ACH 


GTEP 






r 

I FPA’'/ = 

1 

» no 

NOT 

PR I N 

T 

RAW 

Si 

JYC and 


SOU 

A 

R 

ES 

r 

IFA\/-^ 

1 

%no' 

NOT 

PRIM 

T 

AV^R 

Ai 

GFiS 






C 

IFRFSD= 

1 

» no 

NO T 

OR J M 

T 

Rfg I 

O' 

UAL SUM 

S 

so 

U 

A 

RES 

r 

I FCOFN = 

1 

» no 

NOT 

PR I M 

T- 

PART 

I . 

AL CO^“E 

F 

I C I 

F 

N 

TS 

r 

I fprfd= 

1 

.DO 

NOT 

TALC 


PRFD 

I ' 

CTFD 7A 

L 

UFS 




r 

I fcnst= 

1 

.DO 

NOT 

HA VF 


CONS 

T 

TERM I 

N 

EO 

u 

A 

T I ON 

r 

I FFAT A= 

1 

.no 

NOT 

PR IN 

T 

I NPi j 

T 

DATA 






r 

I FPWR= 1 

« 

no 

NOT 

CONS 

TO 

)£=:r V 

Al 

R I ARLES 


TO 

2 

NO POU/FR 

C. 

IEPLOT= 

1 

. no NOT 

PLP 

T 

data 









I japf=6LTAPfQQ 
I N ( 1 ) = 1 OHl POE HGMH 
n I v=SHAMO 

PR0B=4H PRR 
YAX=6HY CALC 
XAX=6HY DATA 
I NVAOM I = I NfVA P- 1 
Kjn 1 M = n 
VAP=n, 

K =n 
FLFVP^L = n . 

NO- NT=n 

njom I M= n 
M0MAX = '^ 


Col, no. 

66 

68 
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APPENDIX C - Continued 


RMnr)ATA = Mn"^A TA 
5 1 = • 

SP = n , 


= • 

I F ( I P-p»„fR ) pnn , 4 1 n , 4 1 
410 NOVAR=?^ ( INVAR- 1 ) +1 
NJVPl =NOVAR+l 
NIX=1 

I F (NOVAR-FO ) 110 , 110,1301 

4 1 F movAR= I NVAR 
NVPl =NOVAR+l 
NT X= 1 

P IF(INVAR-FO) 110,110,1300 
110 no IPO 1=1 ,NV^1 
1 or no IPO j= 1 , NV'^’ 1 
1 PO VFCTOR ( I , J ) =0 
140 IF ( I '^WT 1 9 0"^ , FOO , 1 *^0 

WR T Tt^- f 3 , 0^*=: ) 

GO TO ^10 


1 FO 

(!•••• 


4 1 4 


4 1 6 
4 1 8 
1 81 
1 Po 
poo 
21 o 
220 
1 PO 
1 70 

r 

( • • • • 


P o 


4PP 
4P4 
1 8P 
F 30 
FFO 
c=^FO 
F40 
F 1 o 
C 
C 


no 1 ■PO N= 1 , NOnATA 

READ INPUT DATA (W/0 w/FIGHTS), FORMAT lo 
READCF, 10 ) (DATA (L) ,L=1 » INVAR) 

1 F ( I ) POO , 4 I 4 ,4 1 P 

data (NOVAR ) =DATA ( INVAR ) 

no 4 16 L = 1 , I NVARM I 

K<= T NVARM T +L 

nATA (!<K )=DATA 

DO 1 PI L = 1 ,NO\/AR 

OATA o ( N ,L ) =DATA ( l ) 

no IPO 1=1 , nOvaR 

VICTOR ( I 1 NOVAP+ 1 ) =VFoTOR ( I ^ NOVAR+ 1 ) +DATA ( I ) 
no PPO J= I , NOVAR 

VFCTOR ( I ^ J ) = VFCTOR ( ] ,J)+DATA(! )*DATA(J) 

COMT I NUP 

VFCTOR ( NVP 1 , NVPl )= VECTOR (NVP 1 ^NVPl ) 4-1 *0 
GO TO P 6 P 

calculation sump variable Wt^IGHTP 

DO F 1 O ^^= ] , nooat A 

READ INPUT OA^A (W, h/i=1GHTF), FORMAT 10 
PPAO ( F , 1 0 ) ( OATA ( L ) , L.= 1 , I NVAR ) , WHT 

I F ( I i^L^'WR ) QO O , 4 p^ , 4 p^ 

OATA (N0VAR)=0ATA ( INVAR) 

I NVARM I = T NVAR- 1 
DO 4P2 L= 1 t I NVARM I 
KK= I MV ARM I +L 

r)ATA(’<K)=OATA(L>*-^P 
no 1 8P L= 1 f NOVAR 
nATAP(N,L)=OATA (L) 

DO F40 1 = 1 , NOvAR 

V- CTOR ( I ^ NOV art 1 ) =VP^TOR ( I ,NOVAR+ 1 ) 4 -DATA ( I ) *WHT 
no B40 j= j , NOVAR 

VrCTORC I f J)=VFCTOP( I , J)+OATA ( I )*DATA (J)*WHT 
VFC TOR ( NVP! ♦ NVP 1 ) =VF rTOR ( NVP 1 « NVPl ) +WHT 

CCMPLFTFO sums ^-CUAR'^'5 ANS) CROSS PROOUCTS. THFSE ARE IN 

STORAGE IN LOCAT I 0N» VECTOR ( ; « J ) . '^HESt: WILL BE PRINTED OUT ON 
TAPF 6 UNDER CONTROL OF STATEMENT 100 
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APPENDIX C - Continued 


565 COMTIMUE 
NOVM I =nJOVAR- 1 

566 NO\/PL = NOVAR+l 

567 WRITE (6»90) NOPROB ♦ NOOA TA < NCVAR , VECTOR ( NOVPL ♦ 

INOVPL ) »EFIN,ffOUT 

IF ( IFDATA ) 900 ,40n,STn 

4^^ WRITP^ (6»‘^) 

no 1 R3 N=lfNOnATA 
r • • • • WR I TE I NPUT data « FORMAT 1 2 

183 WRITE(6*12) ( O AT A S ( N , L ) » L= 1 * I NVARM I ) ♦ DAT A<=; ( N » NOV AR ) 

WRITF(6^ 1 1 ) 

570 IF (IF RAW) 900,580,6^0 
580 WRITE (6,15) 

590 WRITE (6,20) ( I , VEC TOR ( I f NOVPL ) « ! = I , NOVM I ) 

600 write (6,25) VEC TOR ( \»0 VAR , NO VPL ) 

610 WRITE (6,30) 

620 WRITE (6,35) ( ( I , J , VFC TOR ( I , J ) * J = I , NOVM I ) , I = 1 , NO VM I ) 

630 WRITE (6,40) ( I , V- CTOP ( I « NO VAR ) , T = 1 , NO VM I ) 

640 WRITE (6^45) VECTOR ( \lO VAR , NOVAR ) 

GO TO 650 

C CALCULATION OF RfsIOuAL FUME OF SOUAREf ANQ CROSS PRODUCTS 

650 IF (IFCNST) 900^651, -^35 

651 IF ( VECTOR ( NOVPL , NOV^L ) ) 652, 652, 65^ 

650 WRIT^^ (6,6"^4 ) 

GO TO olO 

655 DO 660 1=1 ,NOVAR 

6~70 no 660 J=T, NOVAR 

660 VECTOR ( I , J)=VECTOR( I ,J)-( VECTOR! I , NO VPL ) *VEC T OR ( J , NOVPL) 

1 /VECTOR (NOVPL ,NOV^L ) ) 

680 DO 690 1=1, NOVAR 

690 AVF ( I ) = VFCTOR ( I ,N0VPl 5 /VECTOR (NOVPL , NOVPL ) 

70 0 IF CTFAVF) 90 n,Tio,TT‘= 

710 WR I TF ( 6 , *=^0 ) 

720 WRITE (6^20) ( T , A VF ( T ) , I = 1 , NO VM I ) 

T30 WRITE (6,25) \/^ (NOV.aR ) 

7 35 IF (IFREFD) 90O,-74O,-7Rn 
740 WRITE (6,55) 

750 WRITE ( 6 , 35 )(( I , J , VFr TOR C I , J ), J= I , NOVM I ) , I = 1 , NOVM I ) 

760 WRITE (6,4^^) ( I , VFC. TnR ( I , NO V AR ) , I = 1 , NOVM I ) 

779 WRITE (6,45) VECTOR ( MO VAR , NO VAP ) 

789 moftep=-i 

78 1 ASSIGN 1 -^PO TO NLIMREo 

782 0FFR = VFCT0R ( NOVPL , NOv‘=^L ) ~ 1 *0 

790 DO 890 1=1, NOVAR 

791 IF (VECTOR! I , I ) ) 792,-/94,810 

79? WRITE (6,793) I 

GO TO 919 

793 FORMAT (31H ERROR REc^IOUAL SQUARE VARIABLE I4,31H IS NFG A T I VE , PROB 
ILEM T*^RM I NAT'^'O ) 

794 WRIT':^ (6,795) I 

796 sigma ( I ) = 1 

797 GO TO ROO 

795 format (IHOlOH VARlAnLE I , 1 3H IS CONSTANT) 

810 S I GMA ( I ) = SORT ( VECTOR r I , I ) ) 

800 VECTOR ( I , I ) = 1 .0 


820 

no 

870 

1=1, NOVM I 

840 

IPl 

= 1 + 1 


84 1 

no 

870 

j= I PI , NOVAR 


850 VECTOR ( I , J>=VECTOR( I ,J)/(SIGMA( I )*SIGMA( J) ) 
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APPENDIX C - Continued 


8 ?^ V8rTOR( j, I ) i I ^J) 

86 '^ IF 90^ , 87'^ , ^ 

87^ WPIT*^ (6f60) 

874 N!0\/'^P = N0\./M I - 1 

pvc; no r = i,\tov^iF 

8 pn T PI = 1 + 1 

888 WQTT*^^ ( 6 , 88 ) ( I , J » VFr tor ( I , J ) » J= I P 1 » NOV’V’ I ) 
aqn v«/R I T«^ { 6 t 4'"' ) ( T , V‘='^TOP ( I ^ NOVAP ) f I = I ^ NOV^ 1 ) 

1 '^r.r moctfp = N^8T8P+ 1 

3 0ni IF C VICTOR ( NOVAR ,NO\/aP ) ) 1 n ^ 8 , 1 nop , i ni ^ 

I nop ^JFTP^^ 1 =N98TFP~ 1 

WRTT<^ (6»1'^04) NnTP^fll 

GO I 781 

1 C 1 n SI GY= S I GMA (\'0 \/AR)*SQpT (VFPTCR(MD'^AP* NO VAR ) /OFFR ) 
1011 IF ( 1 • 8~vFCT0R (NOVAP ,N0VAR ) ) 1 0 1 3 ^ 1 0 1 2 i 1 0 1 ? 

1 n 1 P Q = 80PT ( 1 • n-V= TOR ( N0\/AR ^ NOVAP ) ) 

GO TO 1017 
1 n 1 n 9 = 0 , 

] \ c=: = 

17 fO^oR) 1 o 1 7 , 1 ni -7 , ^ non 

|0|v tn'Rirr (6^1010) NJ08---P 

GO TO 1 881 
1 Opo I N = o , n 
108 ^ VMAV=^*o 
1 088 MO I N = 0 

1 o^n no I = l,MOVf--^I 

104 1 IF ( VFGTop ( I * I ) ) : 040 , 1 080 , 1 060 

1 04P WRI-TP- (8,1044) I ,NOCx<^^ 

I 04 c^ GO TO 1881 

106^ IF ( V'^CTOP ( I » I } -'^OD 1 0'="O , 1 OBO , 1 ^80 

1080 VAP = VICTOR! I < NOVAP ) * v" CTOP ( NOVAP i I )/VECTOP( I 4 I ) 

I OQn it=- (\/AP) 1 1 4 1 ^7^ , 1 1 1 o 

lion NO I N = NO I M + 1 

1 1 po I MO- y ( NO I N ) = I 

1 130 COFN(NOIN)=VFCTOP ( I , NO VAR ) *G I G^i A C NO VAP ) /S I GM A ( I ) 
1140 SI O ( NO I N ) = ( S I T OMA ( I ) ) ^fqrt C VFCT OR ( I 4 T ) 1 

118^ IF (WIN^ 1 1 6" 4 1 1 4 

904 \nQiT^ (6 40 "' 6 ) 
r:o TO 010 

1 1 80 \/v I N = \/ aP 

1 1 8^ NOM T N= I 
IlQn GO T-o ^nc;n 

116^ 17 (VAP~VMIN) 1 , 1 oc;o , 1 I -rn 

1 1 1 n I cr ( v/A P-\/MA X ) I'^FO, 1^^04101^ 

1 p 1 O \/M Ay = \/ AP 

1 ppo NO'^''Ay=I 

1 08^ CONT 1 Ni J7 

1 ?80 IF (N'OIN) 90 o , 1 p^^ ^ j ^4 .'_5 
QO 8 WPIT7 (64 0'^‘7) 

GO TO 910 

I P4^ WRIT7 ( 6 , 68 ) ''-TGV 
1 P 6 O GO TO 1 800 

1 P48 17 (irCNGT) , 1 OO"'' , 1 P4 6 

1 P46 GMCt = '^ • ^ 

1P4^ GO lOO'^ 

1 p70 gNST = A 'x/7 ( NOVAP ) 

1?7^ no 1880 1 = 1 , N'^ T N 

1 pQo j= I Noi^y ( I ) 

I pq-^ FNFT=ONJQt_ ( nOrr^4 ( r ^ ( J ) ) 
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APPENDIX C — Continued 


1 300 IF (IFSTF^P) , 1 -11 , 1 

1310 IF (NOFNT) 1 31 1 O ■^l 1 , 1 *^1 3 

1311 WRITF (6^91 ) NIOST'^P,!^ 

1^1? GO TO 1 31 

131^ WPITF (6,03) NOF-t^o,.^ 

1314 WRTTP^ (6,70) FLF V*^L , I GY , P , CNST , 

1 { I NOFX ( J ) , CO^N ( J ) , G I r.-'l-O ( J ) , J= 1 , NOT N ) 

131F GO TO ( 1 316, 1 fro ) ,NIy 
1316 GO TO NU^^RFR, (1330,1^^80) 

1 330 FLFVRL = VM I FP/vrr T-QP ( NOVAQ , NOVAP ) 

1 330 IF (FFOUT + FL---VFL) 1 3'^C , ! gro , 1 34 0 
1340 k=NOM1N 
1345 NO'^NT = 0 

GO TO 1391 

1 350 DFN0M = VECTOR ( NOV AP , No V AR ) - V^: AX 
IF (OrNQV) 1 351 , 1 , 1 

1 3Fi NIX=P 

GO “O 1 3TO 

1 3FP FLFVFL = V'^AX*Of^f^P/0‘^NOP/ 

136"^ IF I N-^L- VRL ) 1 3 , 1 36 1 , 1 3RO 

1361 IF (REIN) 1380,1-500,13-0 

1 3 TO l^ = NOVAX 
] OQO N0'^NT = K 

lORl IP (K) 1 ^OP , 1 OQP , 1 4 O 0 , 

1393 WRITF (6,139F) nogt^o 
1 0 Q 4 GO -o Q 1 ^ 

14^0 no 1410 i = i,no\/ar 

14PO IF (I-K) 1 43^^ , 1 4 1 ^ , 1 4 '’o 

1430 no 1440 J=l, NOVAP 

1450 IR (J->^) 1460,144^,1460 

1 4 60 VFrTOR ( I , J ) =V~OTOR ( T , J ) - ( V'^^f OR (!,!<') ^-V^‘CTOR (< , J ) / VECTOR ( K ,< ) ) 

1 440 r OK'T I Nt i'=^ 

14 1'^ FONT I NWF 

147^ no 145^ T = l,NOVAf=5 

1490 IF (I-K) 1 F-'" , 1 4R''’ , 1 P'^O 

Y 500 V'^^CTOR ( I , K ) =-V^C*^OR ( I , K ) /\/R<^ TOR ( l<T , i<r ) 

1 4 RO rONTTNN'^^ 

IFIO no IFP'^ J=!,NnvAP 

1530 IF (J~l<) 1 54n , 1 FP'- , t -:4o 

1 540 VECTOR ( K , J ) = \/FCTOR ( '< , J ) /\/^^ T^r { i< , i< ) 

1 *=^30 rON“ I NUF 

1 FFn VO'OTOR ( K ,K ) = 1 , '^ / V^^^ T pR ( '< , ) 

1 «=^6n 00 TO ( 1 ^ 00 , 1 1 J J V 

1561 WR 1 TF ( 6 , 1 0 P 4 ) nOGTc^P 
R=1 .no 

s T OY=:0 • o 
I FGTFP=0 
GO TO 1 n 1 F 

1 PBO WRITF ( 6 , 7--) nog Tern 
1 3P1 IF (TFc;yppj O n o , ] F O 0 , ] 7 O 

1570 ASSIGN 1 5R0 TO N’JVFcrp 
1 F-7 1 OO T^ 1 1 

1 5R0 FONT I Ml JF 

WRITF (6,15R6) (L , VFrTOR (L , L ) ,L= 1 ,NOVf^ I ) 

1 «=^R 1 IF (TFPPi^n) ooo , 1 c^R p , o 1 n 
1 FRP COMTTNW*^" 

1553 WRITF (6,85) 

1 FP6 format (pAHO OT agonal flfv.fNTS //pOiH VAR .NO* VALNF// 

1 ( 1 H I 7,«=^16.6 ) ) 
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APPENDIX C - Continued 


1 no N=1 ♦NJOnATA 

1^10 YPRPn=rN9T 

] no i = i#NjoiM 

1 K= I NHPX ( T ) 

I 6. -an YPRpn= vpRpn+ro'^N ( i ) -J^nATA p ( n ^ i<) 

Y ( M ) =YPPf^n 
X (N» > =OAT/iS { N #MOVAP ) 

1 npV = f^ATAP (N # NOVAP )-YnRP^n 

PCTDPV = 0'^V*1 '^O./OATA'^ (N^NOVAP) 

Sl=Sl+AR<^ COATAP(NJ f NJOx/AP ) ) 

< 5 P = QP + ARP ( YPPPn ) 

ST = <r,T + /^pc; ( n^V ) 

SA = PA + APP (PrTnc=*V) 

1 WR ITP ( 6f PO ) PATAP C\i ,is)OVAR ) , vnRPQ , PPV » PCTPEV 

P] =S 1 /RNPHATA 
pp = pp/RisjnnAT A 

S-a=pT/RNJOn AT A 
SA = P4/PNtnnATA 
WRTT'^(6^R1 ) f , ^4 

Tf^fT^PLO“r) Q'^n,iAS'^P,nin 
1 eP^ RNjronP( 1 0, IN (^ > ) m V^'^PORiNOPPOR 

1RP6 PORMAT (ARfAA,!'^) 

ENCOnE(6» 1 6P7, YAXI^. ) VAX 
PMrnnP- (6. 1637,XAXI^ > XAV 
1637 PORmaT (A6) 

CALL DD I PLT (1 , I N ♦ NOnA T A . X , Y i 0 » 0 ^ 0 , 0 » 1 ♦ XAX I S , 1 ♦ YAX I 1 ♦ I TAPE ) 
910 CONTINUE 
1 Pi no f n 1 o^' 

18^0 t/(fpTT'=^ (6*90 ) tkjvaR 
U/PfTR (6*94) N9VAP 
no TO QIC 

format 0.‘T,OTO, 1 H lOJO) 

9 format (33X16HT0TAL |NP'JT OA-tA // ) 

C*...P0RMAT 10 P9R Pf^AOlNO I NP'JT DATA 


]n fOPM'^T (nFlP,R) 

) ) prnRM AT ( 1 H^/ ) 

r..*«POR'VAT IP ROP WPITINO INPUT OATA 
) p format ( ^p- 1 p , f ) 

IE PORmaT CIH 49 M 

PO PORMAT (1H IIH SMV X( I?*3H) 

18H «^UM XC IP*OH) =t^iP.4*RH SUM 


SUM OF VARIABLES//) 
= f^iP.4*RH SUM X( IP*3H) =P12.4* 
X( IP*3H) =P^1P.4 ) 


PF format ( 1 *7H si im V =P1P*4) 


on format ( 1 HO 70 H 

IMP CROSS DROnUCTS// 








) 


3E p^ORmat 
1 

P 

40 format 
1 

p 

4F format 
Fn format 


( 1 H 
6 H 
6 H 
( 1 H 
6 H 
6 H 
( 1 H 


7H 


7H 


XiJPn^H) VS X(IP,3H) = P15 

X(IP*th) vs X(I2*3H) = 1 E • 6 , 

X(I2*th) vs X(I2*3H) = PIE. 6 


2 1 M 


( 1 H^ 6 PH 


1 vapia°lps// 

Ft^ forma T ( I H07TH 


X ( 12 ^ 1 ) 
X ( I P * 1 PH ) VO 
X ( I P * 1 PH ) VS 
Y \/S 

) 


VS 

Y 


Y =^1E. 

= >^1 E.6* 

= '=-l ‘=^.6 ) 

= P 1"^.6 ) 


) 


6 » 


average 

R<^SIOUAL FUME 


value op 
OP SO'M 


IRFC, aNO CROCE PROni lOT-^// ) 

60 format ( 1 H069H simple CORRELATI 

1 ON CC^PP I C I PNTS// ) 

6E format (P^^-'HO STANnARO =^9909 OP Y = PI P.6 ) 
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APPENDIX C - Concluded 


70 FORMAT (IIH F LFVF(_ Fl?.4 STANOARO FRROR OF Y F12.4/ 

1 35H MULTIPLF CORRFLATION CCEFFICIFN'T Fl2,F/ 12H 

2 CONSTANT F13*E/^6H VARIABLE COEFFICIENT STD ERR 

30R OF COFF // ( 1 F,H X-T3*F1F.E, F 1 )) 

7F c-q^mAT ( 1 oh COMPLFT^n IF,POH STE^^F OF R‘^GRpr I qM j 

BO POR<^AT (lAX^PH '^1P«F,3H Fip*B»?M Fi7.f=:,3H 

81 FORMAT (IH^flXlPH A\/fRACFF,2H FIP.F^-^h F12.E,2H Fip.E»PH 

1 FI P.'=^ ) 

8E FORMAT (1HP43H PREmcT'^D VS ACT'JAL RF^ultc; /73H 

1 actual PRFniCTFO OrvIATIOM PCT* OF 

PVm //) 

90 format ( 2?H 1 STFPW I SF RfgRFSSION //IPH PROpLFM uq 110 // 1 3H NO 

IDATA = IF //1 8H NO OF VARIABLES = 110 //30H WEIGHTEn OFGRff.c; qf fR 

2FEDOM = F12.2 //?BH LFVFL TO ENTftr vARIARLf = F10.3 //20H f lfvE 
3L TO RfmovE VARIARLF = F9,-^ /// ) 

91 FORMAT (9HOSTFP NO* If / 1 9H VARIABLE REMC'/EO 18) 

92 FORMAT (9H0STEP NO* If /20H VARIABLE FNTFRING 18) 

93 format (4BH1 INVAR SHOULO NOT BE GRfATFR THAN BO* IT IS IB) 

94 format (4BHINOVAR SHOULD NOT BE GRFATFR THAN BO* IT is IB) 

6R4 format ZFRO N' 'mB^R OF DATA* SO LONG*) 

90B format C42H FPROR IN CONTROL CARD^ PRORL‘=‘M TFRMINATfo) 

906 FORMAT (2BH frrcR » V^IN PLUS* SOLONG ) 

907 format (25H '^RROR^NOlM Ni 'S « SOLONG ) 

1004 FORMAT (1H037HY SOUAPR NON-POS I T I VF ^ TE RM I N A T F STEP I B) 

1019 format (IH'^PRH NO MOOT OEGRFES PRi:r^r^O[V' ST'=^'^ I B ) 

I 044 FORMAT ( 1 HO l OH SOUARf X- I B » 1 7H NEGATIV^ * SOLONG I B ^ 6H '^TFRS) 

1 39B format ( 1 2H K = 0, I6, th SOLONG) 

F\in 
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TABLE I.- MATERIAL-PROPERTY AND ABLATION -PERFORMANCE DATA 


Graphite 

Bulk 

density, 

Pp, 

g/cm3 

’ Air dis- 
placement 
density, 

^air» 

g/cm^ 

Theoretical 

density, 

g/cm^ 

Total 

pore, 

volume, 

Vp. 

cm3/g 

Most prob- 
able pore 
radius, 
Rp, 

Mm 

Maximum 
grain size 

mm 

Pore 

volume 

Vp, 

cmVg 

Pore 

area 

Ap, 

mVg 

Surface 

area, 

s, 

m2/g 

Average 

pore 

radius, 

Ta» 

A 

Most probable 
micropore 
radius, 

^P» 

A 

Mean 

pore 

radius 

A 

AHDG 

1.877 

2.041 

2.257 

0.066 

^1.77 

0.84 

0.0033 

1.2 

1.4 

100 

15 

49.4 

ME 11 

1.627 

2.114 

2.255 

.157 

2.36 

.076 

.0025 

1.0 

1.2 

108 

15 

43.6 

ME 14 

1.735 

2.067 

2.256 

.091 

1.74 

.076 

.0027 

.9 

1.0 

120 

14 

53.7 

ME 15 

1.739 

2.077 

2.245 

.091 

2.01 

.18 

.0029 

1.2 

1.3 

100 

15 

45.7 

MEOB 

1.592 

2.15S 

2.255 

.149 

2.42 

.076 

.0026 

.9 

1.1 

109 

15 

49.2 

H-205 

1.783 

2.095 

2.251 

.088 

6.55 

.41 

.0021 

.8 

.8 

120 

14 

52.3 

H-205-85 

1.835 

2,053 

2.248 

.069 

4.74 

,15 

.0023 

.8 

.8 

124 

15 

56.0 

MHLM 

1.777 

2.089 

2.259 

.095 

15.2 

.84 

.0021 

.7 

.8 

127 

15 

53.0 

• MHLM -8 5 

1.857 

2.043 

2.259 

.082 

^2.13 

.84 

.0025 

.9 

1.0 

111 

15 

50.1 

2BE 

1.514 

2.159 

2.257 

.205 

2,13 

.15 

.0029 

1.1 

1.2 

112 

15 

50.6 

2D8D 

1.480 

2.074 

2.230 

,196 

1.77 

.18 

.0042 

1.7 

1.8 

102 

15 

47.6 

2D9B 

1.689 

1.974 

2.252 

.097 

2.13 

.18 

.0032 

1.3 

1.4 

101 

15 

46.6 

W119 

1.739 

2.078 

2.257 

.095 

1.16 

.25 

.0054 

1.9 

2.1 

118 

15 

54.0 

L-56 

1.607 

2.070 

2.232 

.142 

3.55 

.15 

.0024 

1.0 

1.1 

119 

9 

50.3 

L-56-GP 

1.603 

2.075 

2.234 

.144 

3.22 

,15 

.0024 

.9 

1.0 

125 

12 

51.5 

P-3W 

1.661 

2.071 

2.257 

,137 

1.84 

.15 

.0049 

1.5 

1,7 

138 

15 

60.5 

P-3W-GP 

1.640 

2.058 

2.258 

.136 

1.87 

.15 

.0043 

1.3 

1.5 

134 

14 

61.7 

E-24 

1.566 

2.100 

2.232 

.157 

3.88 

.13 

.0027 

1.1 

1.3 

76 

16 

42.5 

3499 

1.680 

2.110 

2.259 

.135 

3.22 

,076 

.0040 

1.4 

1.6 

109 

17 

52,9 

3499 -S 

1.590 

2.183 

2.258 

.177 

4.09 

.076 

.0056 

1.7 

2.0 

133 

15 

58.6 

39 -RL 

1.631 

2.189 

2.258 

.159 

2.84 

.076 

.0051 

1.8 

2.0 

115 

15 

53.2 

4007 ! 

1.661 

2.158 

2.259 

.150 

5.35 

.20 

.0063 

1.9 

2.1 

' 139 

15 

60.7 

8827 

1.764 

2.120 

2.259 

.107 

2.53 

.076 

.0038 

1.1 

1.2 

143 

15 

63.7 

9-RL 

1.667 

2.212 

2.259 

.142 

3.81 

.076 

.0062 

1.9 

2.1 

1 134 

15 

60.9 

9050 

1.777 

2.090 

2.258 

.097 

3.44 

.076 

.0024 

.9 

1.0 

104 

17 

51.3 

LI 

1.539 

2.105 

2.255 

! .192 

2.73 

.15 

.0031 

1.2 

1.3 

114 

16 

52.4 

L31 

1.606 

2.075 

2.231 

.136 

4.63 

.15 

.0016 

1.1 

1.0 

94 

8 

41.9 

331 

1.705 

2.084 

2.257 

.106 

1.19 

.076 

.0035 

1.0 

1.2 

134 

15 

57.7 

AGSX 

1.645 

2.173 

2.259 

.152 

1.69 

.41 

.0033 

1.2 

1.4 

113 

15 

50.6 

ATJ 

1.717 

2.062 

2.257 

.112 

3.55 

.15 

.0020 

.7 

.6 

91 

13 

43.8 

ATJ-GP 

1.709 

2.070 

2.258 

,100 

2.96 

.15 

.0030 

.9 

1.1 

118 

14 

52.3 

ATJS 

1.837 

2.052 

2.257 

.076 

^1.06 

.15 

.0022 

1.0 

1.0 

101 

15 

44.9 

ATJS-GP 

1.825 

2.026 

2.257 

.077 

^1.11 

.15 

.0035 

1.1 

1.3 

123 

15 

57.8 

ATL 

1.791 

2.077 

2.258 

.089 

^1.06 

.76 

.0023 

1.0 

1.0 

98 

15 

46.3 

ATL-GP 

1.776 

2.073 

2.258 

.091 

^1.63 

.76 

,0028 

1.0 

1.1 

124 

13 

54.0 

CDA 

1.654 

2.195 

2.257 

.150 

4,34 

.15 

.0024 

.9 

1.0 

121 

15 

52.3 

CDG 

1.503 

2.151 

2.255 

.214 

6.45 

.41 

.0029 

.9 

1.0 

124 

15 

56.3 

CDG-GP 

1.491 

2.284 

2.257 

.239 

5.35 

.41 

.0061 

1.4 

1.6 

156 

18 

76.7 

CMB 

1.810 

1.978 

2,257 

.044 

.43 

.076 

.0026 

1.1 

1.2 

107 

12 

47.0 

PGR 

1.688 

2.096 

2.259 

.122 

^1.33 

.76 

.0030 

.9 

1.0 

143 

14 

58.5 

AGOT 

1.736 

2.182 

— 

.115 

2.13 

.76 










AXF-5Q 

1.803 

2.208 

— 

.095 

.44 

.025 











Graph -i -tite "A" 

1.897 

2.165 

— 

.050 

^.54 

.84 








___ 


Graph -i-tite ”G" (0.840) 

1.908 

2.050 

— 

.060 

^.54 

.84 










Graph -i -tite "G” (0.203) 

1.873 

2.015 




,203 I 


— 


— - 

... 

--- 


^Best estimate. No definite most probable size. 

^Ncne detected. 

c Specimen lost during run. Length estimated from film record. 
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TABLE I.- MATERIAL-PROPERTY AND ABLATION -PERFORMANCE DATA 


Concluded 


Graphite 

Total 

void 

porosity, 

% 

Air -open 
void 

porosity, 

®air» 

% 

Mercury - 
open void 
porosity , 

Y 

Ash 

content, 

A, 

ppm 

Interlayer 

spacii^g, 

dc, 

A 

Thermal 

conductivity, 

K 

W/cm-K 

Specimen 

length 

change, 

Al, 

cm 

Surface 

tempera- 

ture, 

T, 

K 

Recession 

rate, 

cm/sec 

Surface 

texture, 

T 

Degree of 
deformation, 

r 

AHDG 

16.8 

8.0 

12.4 

1132 

3.361 

1.444 

0.3686 

2178 

0.01432 

2 

2 

ME 11 

27.9 

23.0 

25.5 

417 

3.364 

.641 

.5072 

2189 

.01726 

1 

1 

ME 14 

23.1 

16.1 

15.8 

254 

3.363 

.694 

.4298 

2161 

.01586 

1 

1 

ME 15 

22.6 

16.3 

15.8 

123 

3.379 

.465 

.4486 

2161 

.01715 

1 

1 

ME 18 

29.4 

26.1 

23.7 

362 

3.364 

.523 

.5222 

2183 

.01802 

1 

1 

H-205 

20.8 

14.9 

15.7 

1519 

3.370 

1.226 

.4356 

2206 

.01763 

1 

2 

H-205-85 

18.4 

10.6 

12.7 

4150 

3.375 

1.119 

,3881 

2178 

.01707 

1 

2 

MHLM 

21.4 

14.9 

16.9 

742 

3.358 

1.454 

.6495 

2233 

.02576 

3 

2 

MHLM-85 

17.8 

9.1 

15.2 

692 

3.359 

1.573 

.5197 

2239 

.02310 

3 

2 

2BE 

32.9 

29.9 

31.0 

1119 

3.361 

.683 

.5626 

2228 

.02001 

1 

2 

2DBD 

33.6 

28.6 

29.0 

2491 

3.403 

.184 

.7866 

2283 

.02003 

1 

2 

2D9B 

25.0 

14.4 

16.4 

3197 

3.369 

.423 

.5596 

2261 

.01782 

1 

1 

W119 

23,0 

16.3 

16.5 

4226 

3.361 

.479 

.5057 

2233 

.01873 

1 

1 

L-56 

28.0 

22.4 

22.8 

(b)- 

3.399 

.446 

.4572 

2100 

.01793 

1 

1 

L-56-GP 

28.2 

22.8 

23.1 

(b) 

3.396 

.460 

.4790 

2144 

.01726 

1 

1 

P-3W 

26.4 

19.8 

22.8 

(b) 

3.361 

.832 

.4295 

2061 

.01734 

1 

1 

P-3W-GP 

27.4 

20.3 

22.3 

3 

3.360 

.836 

.4285 

2044 

.01648 

1 

1 

E-24 

29.8 

25.4 

24.6 

932 

3.399 

.400 

.6172 

2222 

.02227 

1 

3 

3499 

25.6 

20.4 

22.7 

466 

3.359 

1.237 

.4943 

2183 

.01970 

1 

1 

3499 -S 

29.6 

27.2 

28.1 

460 

3.360 

1.076 

.4948 

2122 

.01642 

1 

1 

39 -RL 

27.8 

25.5 

25.9 

(b) 

3.360 

1.162 

.4816 

2128 

.01840 

1 

1 

4007 

26.5 

23.0 

24.9 

380 ; 

3.359 

.942 

,4829 

2172 

.01955 

2 

3 

8827 

21.9 

16.8 

18.9 

419 

3.358 ' 

1.191 

.4097 

2094 

.01733 

1 

1 

9-RL 

26.2 

24.6 

23.7 

(b) 

3.359 

1.191 

.4001 

2083 

.01795 

1 

1 

9050 

21.3 

15.0 

17.2 

391 

3.360 

1.240 

.3947 

2061 

.01454 

1 

1 

LI 

31.8 

26.9 

29.5 

709 

3.364 

.673 

.5278 

2111 

.01912 

1 

1 

L31 

28.0 

22.6 

21.8 

377 

3.400 

.460 

.5415 

2150 

.01767 

1 

2 

331 

24.4 

18.2 

18.1 

1476 

3.362 

.537 

.4917 

2167 

.01836 

1 

1 

AGSX 

27.2 

24.3 

25.0 

580 

3.359 

1.313 

.6284 

2183 

.02668 

3 

3 

ATJ 

23.9 

16.7 

19.2 

1209 

3.361 

.970 

,4559 

2200 

.01718 

1 

1 

ATJ-GP 

24.3 

17.4 

17.1 

25 

3.360 

.908 

.3835 

2139 

.01582 

1 

1 

ATJS 

18.6 

10.5 

14.0 

1625 

3.361 

1.153 

.3985 

2161 

.01533 

1 

1 

ATJS-GP 

19.2 

9.9 

14.1 

(b) 

3.361 

1.175 

.3668 

2061 

.01447 

1 

1 

ATL 

20.7 

13.8 

15.9 

5810 

3.360 

1.037 

.6606 

2300 

.02092 

3 

3 

ATL-GP 

21.3 

14.3 

16.2 

17 

3.360 

1.034 

.7142 

2256 

,03072 

3 

3 

CDA 

26.7 

24.7 

24.8 

709 

3.361 

.812 

.4882 

2189 

.01888 

1 

1 

CDG 

33.4 

30.1 

32.2 

531 

3.364 

.801 

.9385 

2211 

.02826 

2 

3 

CDG-GP 

33.9 

34.7 

35.6 

(b) 

3.362 

.819 

.7435 

2178 

.02235 

2 . 

3 

CMB 

19.8 

8.5 

8.0 

64 

3.362 

.380 

.4178 

2139 

.01554 

1 

1 

PGR 

25.3 

19.5 

20.6 

688 

3.358 

.781 

c 1.0460 

2244 

.04287 

3 

3 

AGOT 

--- 

20.4 

20.0 

*^740 

— 

^1.450 

.5319 

2200 

.02252 

3 

2 

AXF-5Q 

--- 

18.4 

17.1 

dl50 

— 

d<.796 

.4567 

2156 

,02162 

1 

1 

Graph-i-tite ”A” 

— 

12.4 

9.5 

«^<600 

— 

‘*.949 

.3719 

2206 

.01468 

2 

3 

Graph-i-tite ’*G” (0.840) 

--- 

6.9 

11.4 

d<600 

— 

** 1.068 

® .3297 

®2200 

® .01259 

^2 

"2 

Graph-i-tite ”G” (0.203) 

— 

7.1 

— 

370 

— 

‘*1.126 

(f) 

(f) 

(f) 

ei 



^Obtained from manufacturer’s literature, 

^ Also see additional ablation data in table n. 
*See table n for ablation data. 
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TABLE n.- COMPARISON OF Graph-i-tite "G" ABLATION 
PERFORMANCE FOR DIFFERENT GRAIN SIZES 
[Nominal heating rate = 562 W/cm2[j 


Material 

Total recession 
in 30 seconds, 
AZ, cm 

Recession rate, 
§, cm/sec 

Surface 
temperature, 
T, K 

Graph-i-tite "G" 

0.3289 

0.01461 

2050 

(Gm = 0.840 mm) 

.3147 

.01327 

2010 

Graph-i-tite "G" 

0.3358 

0.01458 

1980 

(Gm = 0.203 mm) 

.3366 

.01423 

1970 


.3457 

.01341 

1950 

ATJ 

0.4224 

0.01757 

2120 


TABLE m.- CRITERIA FOR ELIMINATING GRAPHITES 
FROM CONSIDERATION BASED ON F s 1.5 


Maximum grain 
size, Gm, 
mm 

Calculated density below 
which if the actual bulk 
density of a graphite falls, 
it is to be eliminated, 
g/cm^ 

0.84 

2.20 

.76 

2.13 

.41 

1.84 

.25 

1.72 

.203 

1.68 

.20 

1.68 

.18 

1.66 

.15 

1.64 

.13 

1.62 

.076 

1.58 

.025 

1.54 
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TABLE IV.- ABLATION PERFORMANCE PREDICTIONS AT 1 -PERCENT SIGNIFICANCE LEVEL 


Graphite 


cm 

T, 

K 

s, cm/sec 

T 

r 

Predicted 

Measured 

Predicted 

Measured 

Predicted 

Measured 

Predicted 

Measured 

Predicted 

Measured 

Equations (16) to (20) 

AGOT 

0.6952 

0.5319 

2207 

2200 

0.02645 

0.02252 

2.8 

3 

2.7 

2 

AXF-5Q 

.3080 

.4567 

>2123 

2156 

.01328 

.02162 

.8 

1 

.7 

1 

Graph-i-tite ^^A” 

.5715 

.3719 

<2251 





2 

2.4 

3 

Graph-i-tite ^^G'' 

.5607 

.3297 

<2243 

2200 


.01259 

3.0 

2 

2.4 

2 

(0.840 mm) 

j 










Equations (21) to (23) 

^AXF-5Q 

(b) 

0.4567 

(b) 

2156 

0.01604 

0.02162 

(c) 

1 

■H 

1 


^From correlations for graphites with the requirements that r ^ 1.3 (Gm = 0.22 mm) and T ^ 1.5 (criteria in table HI) 


at the 2, 5 -percent significance level. 

^No prediction possible because not available. 
^Not applicable. 
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Graphite test specimen 



Water - cooled holder 



Phenolic- asbestos insulator 


Figure 1.- Schematic drawing showing graphite test specimen, water-cooled holder, and insulator. 



sped lit n test t i «e| 


sec 


Figure 2.- Sample plot of specimen length and surface temperature as a function of time (AHDG graphite). 




Untested 


T = 1 (smooth) 
r = 1 (no gonging) 



T = 2 (sli^tly rouj^) 
r = 3 (severe gouging) 



T = 2 (sli^tly rou^) ■'' = 3 (rough) 

r = 2 (sli^t gouging) r = 2 (slight gouging) 

L-72-116 



Figure 3.- Photograph of representative after-test specimens indicating the quantitative 
assignation of surface texture t and degree of gouging T. 
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(e) Pore volume, equation (5) 
Figure 4.- Continued. 
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